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Using atomistic quantum simulation based on a tight binding model, we investigate the formation 
of energy gap E g of graphene nanomesh (GNM) lattices and the transport characteristics of GNM- 
based electronic devices (single potential barrier structure and p-n junction) taking into account the 
atomic edge disorder of holes. We find that the sensitivity of E g to the lattice symmetry (i.e., the 
lattice orientation and the hole shape) is significantly suppressed in the presence of the disorder. In 
the case of strong disorder, the dependence of E g on the neck width is fitted well with the scaling 
rule observed in experiments [Liang et al, Nano Lett. 10, 2454 (2010)]. Considering the transport 
characteristics of GNM-based structures, we demonstrate that the use of finite GNM sections in 
the devices can efficiently improve their electrical performance (i.e., high ON/OFF current ratio, 
good current saturation and negative differential conductance behaviors). Additionally, if the length 
of GNM sections is appropriately limited, the detrimental effects of disorder on transport can be 
avoided to a large extent. Our study provides a good explanation of the available experimental data 
on GNM energy gap and should be helpful for further investigations of GNM-based devices. 
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I. INTRODUCTION 

In the recent years, thanks to its specific band struc- 
ture and excellent carrier transport properties [H [2], 
graphene has been expected to be an excellent candi- 
date for electronic applications [3H5]. However, the lack 
of energy gap between valence and conduction bands is 
a serious drawback regarding the operation of graphene 
transistors. In particular, the zero bandgap makes it very 
difficult to have a high ON/OFF current ratio and a re- 
ally saturated current at high drain voltage [5J- Fortu- 
nately, a number of techniques to open a bandgap in 
graphene have been proposed, for instance, cutting a 
graphene sheet into nanoribbons [7j , applying an electric 
field perpendicular to a bilayer graphene sheet [S], inter- 
action of graphene with the substrate [10], or intro- 
ducing doping and impurities in graphene [111 112] . This 
has hence stimulated a lot of investigations of various 
graphene nanodevices [T3H22] . 

Besides the finite bandgap graphene-based structures 
mentioned above, experimentalists have recently re- 
ported the fabrication of a new graphene nanostructure 
[23U26] called graphene nanomesh (GNM) in which the 
size of nanoholes and the neck width (shortest distance 
between neighbor nanoholes) can be controlled down 
to the sub-10 nm scale. Various techniques have been 
developed to produce and control such GNM lattices, 
e.g., block copolymer lithography [231 121], nanosphere 
lithography [25] . and nanoimprint lithography [26] . This 
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nanostructuring can open up a finite bandgap in large 
graphene sheets and hence GNM-based FETs can ex- 
hibit a similar ON/OFF current ratio as in graphene 
nanoribbon (GNR) devices [23]. Especially, GNM (2D 
semiconducting film) devices are able to carry electric 
current about 100 times greater than individual GNRs. 
The interest of GNMs is additionally due to the possi- 
bilities of bandgap engineering when varying the neck 
width [53J [55] . The GNM lattices are hence expected to 
enlarge the applications of graphene-based devices, not 
only in electronics. 

In this context, the theoretical analysis of electronic 
structure and transport properties of GNM lattices are 
very timely and desirable to provide a guidance for fur- 
ther investigations. Actually, a number of works [2"71 - f3"?] 
have been carried out in this direction. The studies in 
[2"7H3"5] have focussed on the electronic structure of per- 
fect GNM lattices, especially, on the property of bandgap 
depending on the lattice constant, the hole shape and the 
orientation (i.e., holes along zigzag or armchair direc- 
tions). Some of these works have tried to propose scal- 
ing rules to explain and predict the behavior of bandgap 
when varying the lattice parameters. There are also a few 
studies in which the transport properties of GNM-based 
structures have been investigated [33H37]. In particu- 
lar, the transport properties related to the formation of 
GNM electronic structure were analyzed in [33] , the ther- 
moelectric transport was considered in I34H36] , while the 
design of negative differential conductance GNM-based 
devices was discussed in [37] . However, to our knowl- 
edge, one important issue, the effects of disordered holes 
in GNM lattices, has not been studied yet. Therefore, our 
aim in this work is to clarify this point and we focus on 
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two subjects: (1) the effects of edge disordered nanoholcs 
on the energy gap of GNM lattices and (2) the optimiza- 
tion of the operation of some devices based on GNMs 
taking into account the disorder. The rest of paper is or- 
ganized as follows. In section II, we introduce briefly the 
model and calculation methods which consist basically 
in using the non-equilibrium Green's function (NEGF) 
technique to solve a tight-binding Hamiltonian as in |37j . 
The formation of the energy gap of GNM lattices is dis- 
cussed in section III. A and the transport characteristics 
of single barrier structures and p-n junctions based on 
GNM lattices are presented in section III.B. Finally, the 
conclusion is drawn in section IV. 



II. MODEL AND CALCULATION 

In graphene nanostructures, the electronic properties 
of charges around the intrinsic (midgap) Fermi level can 
be described well by a single orbital tight binding ap- 
proach 27-2{Jj. The corresponding Hamiltonian reads 



h.c 



(1) 



(n,m) 



where U n is the on-site energy/potential energy at the n- 
th site and t = 2.7 eV [3§] denotes the hoping energy be- 
tween nearest neighbor carbon atoms. This Hamiltonian 
has been recently used to investigate the electronic struc- 
ture and the transport characteristics of perfect GNM- 
based devices [37]. In such 2D structures, the extension 
along the transverse (OY) direction is considered through 
Bloch boundary conditions. Using the same approach 
as in [37], the considered lattice is split into elementary 
cells, the operators in equation (1) are then Fourier trans- 
formed along the OY direction and the Hamiltonian is 
rewritten in terms of decoupled (quasi-lD) forms 

H tb = ^2H lD (ky) (2) 

ky 

HlD(ky) = Y,tin(ky)+H n (k y ) + H+(k y ) 
n 

where H n (k y ) is the Hamiltonian of cell {n} and H„ (k y ) 
denotes the coupling of cell {n} to cell {n ± 1} along the 
transport direction OX. 

Using the Hamiltonian (2), the device retarded Green's 
function in the ballistic approximation is computed for a 
given momentum k y as 



G (e, k y 



Hid (k y ) - E (k y ) 



(3) 



where E (k y ) — E^ (k y ) + E^ (k y ) with the self-energies 
El, ft describing the left and right contact-to-device cou- 
pling, respectively. The local density of states and the 
transmission probability needed to evaluate the current 
are calculated as D(E,k y ,f n ) = — lva{Q n n (E 1 k y )} /ir 
and T(e, k y ) — Tr {r^r^C^}, respectively, where 



is the transfer rate at the 
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left (right) contact. The current is then computed using 
the Landauer formula 



J = 



nh 



dkn 



deT(e,k y )[f L (e)-f R (e)} (4) 
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where f L{R) (e) = [l + exp ((e - E FL ( R) ) /k b T)] 1 is the 
left (right) Fermi distribution function with the Fermi 
level Epimy At zero bias, the transport properties of 
the system is characterized by its conductance 
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where the transverse width of graphene sheet Lq is as- 
sumed to be much larger than the device length along 
the transport direction [37] . The integrations over k y are 
performed in the first Brillouin zone. 

The electronic structure and the transport in some 
GNM-based devices have been studied in [2"7H37] but 
most of nanohole lattices are perfectly periodic. In prac- 
tice, the disorder, e.g., the hole shapes may be randomly 
different, is inevitable in the fabrication process of GNM 
lattices [53] ES]- Therefore, using the formalism above, 
we investigate the effects of this disorder on the electronic 
structure and the transport characteristics of GNM struc- 
tures. To model the disorder, the carbon atoms in the 
edge of nanoholes are randomly removed with a uniform 
probability P D [38H40] . We focus on two different GNM 
lattices that consist of holes along the armchair (arm- 
chair holes) [37] and zigzag (zigzag holes) [2S] directions, 
respectively. 



III. RESULTS AND DISCUSSION 

A. Energy gap of GNM lattices 

First, we would like to review the recent studies on 
the electronic properties of GNM lattices. Strong efforts 
have been made on theoretical models and calculations on 
GNMs [27H33] to predict the properties of their electronic 
structure and, especially, the possible scaling rule for 
bandgap E g with respect the shape, size, and architecture 
of holes introduced into the graphene sheet. However, the 
obtained results are very different. In particular, Peder- 
sen et al. |27j considered the lattices of armchair circular 
holes and found a scaling rule E g oc y/N rem /N t ot, where 
N rem is the number of removed carbon atoms in a super- 
cell originally containing N to t atoms. However, this sim- 
ple rule is not applicable for large ^/N rem /N tot when E g 
fluctuates very strongly with the change in such the ratio. 
Moreover, together with the results obtained in [28] . it is 
shown that the behavior of E g is also dependent on the 
hole shape (circular holes in |27j . triangular and rhom- 
bus holes in [IE])- More significantly, the studies on the 
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FIG. 1: Bandgap of GNM lattices as a function of the super-cell lattice constant for different hole shapes, (a) is for zigzag holes 
while (b) is for armchair holes. The super-cells of considered GNM lattices, where white circles indicate the removed carbon 
atoms, are presented in the bottom of two sub-figures. 



lattices of zigzag holes have shown very different results. 
Ouyang et al. [SU] predicted half of GNMs are semimet- 
als and the rest are semiconductors, while it was shown 
in other works that only one-third of considered GNMs 
had a significantly finite bandgap [33 GH] or even all 
GNMs are semiconductors [28] . The authors in ref. (25] 
also reported the scaling rule E g oc 2~Q/ 3 for this kind of 
GNM lattice, where Q is the super-cell lattice constant 
in the unit of a c = 2.46 A. The complex properties of E g 
can be explained by the strong sensitivity of the GNM 
electronic structure to the hole shape and obviously to 
the lattice orientation, i.e., to whether we have zigzag or 
armchair holes. 

Experimentally, Liang et al. |26j reported the bandgap 
scaling rule E g oc 1/W nw (W nw is the neck width), which 
can be explained by the important role on the GNM 
bandgap of the quantum confinement in the multiple 
graphene nanoribbon (GNR) network. A recent theoret- 
ical work [33] reported the same scaling rule for E g when 
studying the lattices of irregular holes but with a quite 
large distribution of E g at each W nw . The latter feature 
indicates that the GNM bandgap is very sensitive to even 
a small change in the hole edges. It is important to note 
that neither metallic GNMs nor the orientation effects 
have been experimentally observed 123H26j . As discussed 
above, the edge disorder of holes in the GNM lattice is 
inevitable [23l [26] . Suggested by the studies [38H40] on 
disordered GNR structures, the edge disorder of holes 
may be an important factor to explain the discrepancies 
between theoretical works introduced above and exper- 
iments. To examine this idea, we model the disordered 
GNM lattices as follows. We assume that when fabricat- 
ing GNM lattices with a given hole shape, a few percents 



of additional atoms at the edge of holes are unintention- 
ally pulled out of the lattice. To model this system, we 
start with a (ordered/original) lattice of perfectly peri- 
odic holes and then randomly remove the edge atoms 
with a uniform probability Pjj, This new (disordered) 
lattice has a mixing of a number of original holes and a 
few percents of holes having different shapes. 

We now use the tight binding Hamiltonian (1) to re- 
compute the bandgap E g of lattices of perfectly peri- 
odic holes to understand its dependence on the lattice 
parameters. In Fig. I, we plot E g as a function of the 
super-cell lattice constant W (i.e., W = Qa c ) for dif- 
ferent hole shapes and for different orientations (zigzag 
holes in Fig. 1(a) and armchair holes in Fig. 1(b)). Gen- 
erally speaking, it is shown that E g of semiconducting 
GNMs inversely decreases with an increase of the super- 
cell lattice constant W . Qualitatively, this is in good 
agreement with the experiments. However, it confirms 
again a big difference between the lattices of zigzag and 
armchair holes. In the case of zigzag holes, the results for 
the hole shape (1) (already studied in [23]) shows that a 
finite bandgap is observed when the index Q is a mul- 
tiple of 3 while all other GNMs are semimetallic. This 
feature has been explained in [29) as a consequence of 
the inter-valley scattering between different Dirac points 
of pristine graphene when the holes are created. How- 
ever, our results in Fig. 1(a) demonstrate this is not the 
case of other hole shapes even with a small change com- 
pared to the hole (1) and, simultaneously, the scaling 
rule Eg cx 2 -< 2/ 3 [29] is not applicable. For instance, the 
lattices with the hole shape (4) behave contrary to the 
previously suggested law, i.e, they are semimetals when 
Q is a multiple of 3 while others are semiconductors. In 
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FIG. 2: (a) Conductance in the unit of Go = e 2 Lc/hW as a function of Fermi energy for different disorder strengths . The 
plot is for the lattices of zigzag holes with the hole shape 1 as shown in Fig. 1(a) and Q — 13 corresponding to the metallic 
GNMs. (b) The energy gap of the GNM lattices as a function of the super-cell lattice constant with disorder effects. The plot 
is for the lattices of both zigzag and armchair holes with the hole shape 1 as shown in Fig. 1(a) and (b), respectively. 



this case, the GNM bandgap decays very rapidly with 
an increase of W. Besides, no semi-metallic GNMs with 
holes (2) and (3) are observed. Considering the lattices 
of armchair holes, we find from the results displayed in 
Fig. 1(b) that all studied GNMs are semiconductors. 
However, E g is also sensitive to the change in the hole 
shape. This sensitivity of E g to the hole shape is in good 
agreement with that reported in [27J I2H] since the ra- 
tio v ' N rern /N to t of the considered lattices is large, i.e, 
> 0.02. Consistently with the previous works, our ob- 
tained results demonstrate again two important points 
for the lattices of perfectly periodic holes: (1) the GNM 
bandgap is very sensitive to the change in the hole shape 
and the lattice orientation; (2) it is hence difficult to de- 
termine theoretically a unique scaling rule to describe the 
bandgap of all different perfect GNM lattices. 

Next, we go to investigate the disorder effects on the 
energy gap of GNM lattices. Here, we choose to present 
the results obtained for the lattices wherein the holes 
numbered (1) in Fig. 1, which are close to the holes 
fabricated experimentally [231 HSJi are used as original 
holes. As explained above, the other holes (e.g., see in 
Fig. 1) probably occur in the considered lattices with 
a few percent fraction. As done in [39l 00] , we use eq. 
(5) to compute the conductance at zero temperature and 
then measure the transport gap around the zero energy 
point, which implies the formation of the bandgap of the 
considered lattices [26]. Note that all transport quanti- 
ties are computed averagely over forty disorder samples. 
In Fig. 2(a), we plot the conductance of the GNM lat- 
tices as a function of energy and with different disorder 
strengths. In this figure, the lattices of zigzag holes with 
Q = 13 are used. Similarly to the metallic GNRs studied 
in [55H1D] . the ordered GNM lattice is a semimetal, i.e., 
no conduction gap is observed for Pjj — in Fig. 2(a). 
When the disorder occurs, the conductance is strongly 
suppressed, especially near the neutrality (zero energy) 
point. Close to the neutrality point, therefore, a deep 



conduction gap develops when increasing the disorder 
strength. The same feature is obtained for the lattices, 
which originate from the semiconducting lattices of per- 
fectly periodic holes, but the enhancement of conduction 
gap is weaker than that in the metallic ones. 

In Fig. 2(b), we display the energy (conduction) gap 
as a function of the super-lattice constant W for dif- 
ferent orientations (zigzag and armchair holes) and dif- 
ferent disorder strengths. In what follows, we use the 
concepts "semiconducting (metallic) GNMs" to identify 
the lattices originating from semiconducting (metallic) 
GNMs of perfectly periodic holes. As discussed above, 
it is shown that while the conduction gap of semicon- 
ducting GNMs is slightly enlarged, the gap of metallic 
ones quickly develops with an increase of the disorder 
strength. Therefore, in the case of strong disorder, no 
metallic GNMs of zigzag holes are detected and the de- 
pendence of energy gap on W tends to a unique scaling 
rule. Additionally, the results presented in Fig. 2(b) 
shows that the orientation effects are suppressed and the 
same scaling rule is obtained for two different orienta- 
tions in disordered lattices, i.e., when Pjj = 8%. More 
important, the energy gap in this case is fitted very well 
with the scaling rule E g — a/W nw [26] with a sa 0.6 
eVnm. The difference between the value of a obtained 
here and a = 0.8 0.95 eVnm in [53] can be explained 
by the fact that our considered holes are smaller and the 
disorder seems to be weaker than in the case studied in 
their work. Besides, our simulations show that even in 
the presence of weak disorder, the orientation effects are 
always suppressed and the scaling rule E g — a/W nw is a 
good approximation for different original holes (not dis- 
played here), though E g (or a) is still sensitive to a sig- 
nificant change in their shape. However, this sensitivity 
of Eg to the hole shape is also suppressed in the presence 
of strong disorder (very large Pd) when all different holes 
contribute equivalently to the electronic properties of the 
sample. 
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These results are very similar to those observed for 
GNRs 7, 38-40] and demonstrate that the edge disorder 
of holes is one of important factors which weaken the 
lattice symmetry effects and make the effects of quantum- 
confinement in the nanoribbon crossing network of GNM 
lattices dominant. It explains well the behavior of energy 
gap observed in experiments |26j . 



B. Transport characteristics of GNM-based 
structures 

In this sub-section, we discuss the use of the GNM 
lattices to optimize the operation of graphene electronic 
devices including the disorder effects. We focus on two 
typical devices based on a single (gate-induced) potential 
barrier structure and a p-n junction. In what follows, the 
zigzag holes are used and all simulations were performed 
at room temperature. 



1. Single potential barrier structures 

The schematic of the simulated single-barrier structure 
is shown in Fig. 3(a), wherein the potential barrier can be 
generated and controlled by a gate electrode as in field- 
effect-transistors [41]. The model used here is not self- 
consistently coupled to Poisson's equation. Though ideal 
and simple, it can provide the basic and most important 
information on the transport behavior and the electri- 
cal characteristics of devices [42] such as the graphene 
FET. One of the important properties of transistors is 
their ability to switch off the current by tuning the gate 
voltage, i.e., by changing the potential barrier. In princi- 
ple, to significantly switch off the current, i.e., to obtain 
a high ON/OFF current ratio, a large bandgap of the 
channel material is essentially required. Therefore, the 
2D monolayer graphene devices have a serious drawback 
inherent in the semi-metallic character of the conduct- 
ing material. Using a GNM lattice with finite bandgap 
[2"31 l2l)] offers a possibility to overcome this limitation 
more conveniently than by using very narrow GNRs {71 
or bilaycr graphene with strong vertical electric field [8]. 

Obviously, the disorder is not only shown to influence 
strongly the band structure, it is also expected to signif- 
icantly affect the transport quantities such as the con- 
ductance and the current [43J [44] . However, to obtain 
a low OFF current, a finite bandgap material is manda- 
tory in the gated part of the channel, but not in two 
source and drain access regions (or leads here). There- 
fore, we propose to use the graphene device as schema- 
tized in Fig. 3(a) with a finite GNM section inserted in 
between pristine graphene leads to limit at the minimum 
the GNM-induced effects of disorder while the good con- 
trollability of the current can still be guaranteed. We plot 
in Fig. 3(b) the current density as a function of the bar- 
rier height Uq for different lengths of GNM section and 



in two cases of ordered and disordered lattices. Here, the 
semiconducting GNM with Q — 15 is used. It is clearly 
shown that the OFF current (ON/OFF ratio) is strongly 
reduced (enhanced) when inserting a finite GNM section 
into the device, i.e., compare the current for Lqnm =/= 
with that for Lqnm — (pristine graphene). Very sig- 
nificantly, in the presence of disorder, the effect of in- 
creasing Lqnm is not only to reduce the OFF current 
and to enhance the ON/OFF ratio, but also to reduce 
the ON current. To have both high ON current and high 
ON/OFF ratio, our results provide the important sug- 
gestion that the optimum device should be determined 
by the competition between the effects of finite bandgap 
and disorder of the GNM section used. The best device 
parameter (Lqnm) is thus dependent on the bandgap 
value and on the disorder strength, i.e., to keep constant 
the ON (OFF) current, the shorter (larger) Lgnm is re- 
quired if the disorder (bandgap) is stronger (smaller). 
For instance, we suggest that Lqnm — 51 -f- 77 nm shall 
be a good choice for the devices studied in Fig. 3(a), i.e, 
the ON/OFF ratio is about a few hundreds while the ON 
current is still high. 

For the devices with a metallic GNM section, the dis- 
order not only results in the suppression of the current 
but also essentially originates the energy gap opening. 
To see the disorder effects in this case, we display in Fig. 
3(b) the current density for different disorder strengths. 
If the holes are perfectly periodic, it is shown that the 
current is smaller than that in normal graphene devices 
(see in Fig. 3(a)). This is nothing but the consequence of 
an enhancement of the channel resistance when inserting 
the nanoholes in the device. However, due to the lack of 
finite bandgap, the ON/OFF current ratio is quite simi- 
lar in both cases. Interestingly, when the disorder occurs 
in the GNM section, though the current is reduced, the 
ON/OFF ratio significantly increases, i.e., it is about 16, 
65 and 310 for P D = %, 4 % and 8 %, respectively. 
When comparing the two cases of semiconducting and 
metallic GNM devices, we observe that the current and 
the ON/OFF ratio tend to the same scale when the dis- 
order is strong enough. This result is very consistent 
with the property of energy gap discussed above. Fi- 
nally, we emphasize again that for both of these cases, 
the ON/OFF ratio may increase but the ON current de- 
creases with an increase of Lqnm in the disordered de- 
vices and hence a finite GNM section shall be mandatory 
to ensure good device operation. 

Regarding the electronic applications, besides the im- 
provement of the ON/OFF current ratio in the graphene 
devices, recent studies [51 |4"5H4"8"] have been done to 
achieve and control a low output conductance in satura- 
tion or even a negative differential conductance (NDC). 
This leads to the improvement of voltage gain, one of im- 
portant figures of merit for analog high frequency appli- 
cations. In gapless monolayer graphene devices the cur- 
rent saturation is relatively weak [6l [41] because of the 
important contribution of band-to-band tunneling [42] . 
In principle, a good current saturation behavior is ob- 
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FIG. 3: (a) Schematic of potential barrier structures studied in this work. A finite GNM section is inserted in the device center, 
(b) and (c) show the current density at V& = 25 mV as a function of the barrier height with the barrier length Lb ~ 61.8 nm 
and Ef = 0.4 eV. In (b), the results are plotted for Q = 15 (W ~ 3.7 nm, semiconducting GNM), for different Lgnm and 
disorder strengths, (c) is for Q = 14 (W ~ 3.45 nm, metallic GNM), and for different disorder strengths. In (b) the case of 
full pristine graphene structure is shown for comparison. 




Bias voltage (V) Bias voltage (V) 



FIG. 4: (a) I- V characteristics of the potential barrier structure with a finite perfect GNM section and for different barrier 
heights, (b) shows the effects of different Lgnm and disorder strengths for the case of Uo — 0.1 eV in (a). The device parameters 
are E F = 0.4 eV, L B ~ 61.8 nm, Lgnm ~ 76.7 nm, and Q — 15. 



served only in transistors made of large bandgap channel course graphene nanomeshes. It is hence important to 
materials such as conventional semiconductors, semicon- analyze the output conductance and the current satura- 
ducting carbon nanotubes, graphene nanoribbons and of tion in our simulated GNM devices and to evaluate the 
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FIG. 5: (a,b,d,e) Local density of states and (c,f) current spectrum at two different biases (14 = 0.2 and 0.4 V) for the device 
studied in Fig. 4(a). The barrier heights are Uo = 0.1 eV for (a,b,c) and 0.6 eV for (d,e,f). The solid line in (a,b,d,e) indicates 
the potential profile. 



impact of GNM disorder, though the GNM section is fi- 
nite and limited to the active device region. In Fig. 4(a), 
we display the I-V characteristics obtained for different 
barrier heights. Here, the semiconducting GNM lattice of 
perfectly periodic holes is used while the barrier length 
is fixed (Lb — 61.8 nm). In the normal graphene de- 
vices, a weak current saturation can be observed and 
explained by the effects of the transmission valleys oc- 
curring around the Dirac (neutrality) points in the left 
contact and the channel [41] [42] . When rising the bias 
voltage, while the extension of the energy window of cur- 
rent transmission [Efr, Efl) tends to increase the cur- 
rent, the two transmission valleys mentioned above when 
approaching the window [Efr, Efl] can reduce the out- 
put conductance. Hence, the output conductance is fully 
dependent on the way these transmission valleys reach 
the window [Efr, Efl]- In particular, the lowest output 
conductance (current saturation or even NDC) should 
be achieved when the Dirac point Eq in the channel 
approaches Efl at the same bias as Efr reaches the 
Dirac point Ef in the left contact. This suggests a way 
to improve such the behavior by using a finite bandgap 
graphene channel to make the transmission valleys wider 
and deeper [41] . It is exactly what we propose in this 
work: using a finite GNM section in the barrier region to 
enlarge the transmission valley around Eq. Indeed, the 
results presented in Fig. 4(a) illustrate clearly this idea, 
i.e., a good current saturation (or even NDC) behavior 
is observed for some potential barriers. 



In particular, the low output conductance can be 
achieved in two regimes: very low and high barriers. To 
explain its origin, we display in Fig. 5 the local density 
of states and the current spectrum as a function of en- 
ergy at two different biases V& = 0.2 V and 0.4 V and 
for two different barrier heights Uq — 0.1 eV and 0.6 
eV. For Uq — 0.1 eV, when rising the bias voltage, Eq 
moves down in energy to come close to and/or even be- 
low Ef. Hence, the transmission valley (wide and deep 
compared to that in normal graphene devices) around 
Eq strongly limits the current spectrum close to E F r 
(probably including the band-to-band tunneling [32]) at 
high bias, i.e., see in Figs. 5(a,b,c). Moreover, besides 
the main energy gap around the neutrality point, some 
mini-gaps emerge at high energy in the band structure 
of GNMs [37]. Here, the first mini-gap in the GNM con- 
duction band generates an additional transmission valley, 
which can reduce the current component close to Efl for 
U$ = 0.1 eV when rising VJ,. Based on these features, the 
current (with a low output conductance) at high bias in 
this case is smaller than for Uq = 0.2 or 0.3 eV. For high 
barriers (e.g., Uo — 0.5 and 0.6 eV), the transmission val- 
ley around E c reach the energy level Efl at high bias, 
which reduces strongly the current component close to 
Efl, while the transmission valleys around E^ and due 
to the first mini-gap in the GNM valence band limit the 
current close to Err, i.e., see in Figs. 5(d,e,f). The 
output conductance is hence small and can be negative 
when rising the bias voltage. For instance, for the barrier 
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Uq = 0.6 cV = 3/2Ep, E c approaches Efl at the same 
bias Vb = 0.4 V as Efr reaches Ef^ and a strong NDC 
behavior is hence achieved. 

To clarify the disorder effects, we plot in Fig. 4(b) 
the I-V characteristics for different Lqnm and disorder 
strengths while the potential barrier is fixed (Uq = 0.1 
cV). First, considering the case of perfectly periodic 
holes, the current is shown to decrease when increas- 
ing Lqnm- This feature is observed because the finite 
bandgap material (GNM) is used not only in the barrier 
region but also in the two contacts when Lqnm > Lb- 
Hence, when increasing Lqnm (> Lb), the effects of en- 
ergy gap in the GNM sections in the two contacts can 
affect and reduce the current, which is consistent with 
the results presented in Fig. 3(a). More important, the 
results for disordered devices in Fig. 4(b) demonstrate 
again that though the low output conductance is always 
observed at high bias, an increase of Lqnm enhances the 
disorder effects, i.e., reduces significantly the current. 

All results presented in this sub-section suggest that 
graphene devices can be optimized by inserting only a 
(appropriately) finite GNM section into the channel to 
avoid the disorder effects while good characteristics such 
as high ON/OFF current ratio and strong current satu- 
ration (or even NDC behavior) are still guaranteed. For 
instance, we propose to use Lqnm — 51 -f- 77 nm for the 
devices considered here. 



2. P-N junctions 

Now, we discuss the use the GNM lattices to improve 
the operation of graphene P — N junctions. The P — N 
junctions are very well-known electronic structures, since 
they can exhibit a non-linear I-V characteristics with a 
strong NDC behavior [2lJ [37] , which offers great poten- 
tial for high frequency applications [IS]- In principle, this 
NDC behavior occurs at high bias when the filled states 
in n-doped side find a reduced number of empty states 
available for tunneling in the p-doped side. The feature 
can be observed essentially due to the large bandgap ma- 
terials used in two doped contacts, e.g., GNRs in [3T] and 
GNMs in [37] . Besides, it was also shown that the current 
peak in this kind of device is strongly dependent on the 
bandgap and the length of the transition region in the 
device center. Therefore, we have proposed in [3T] [37] to 
use graphene hetero-structures with a gapless (or small 
bangap) in the transition region to obtain high current 
peak and peak-to- valley current ratio (PVR). However, 
there is still one issue, the disorder effects, which has not 
been investigated yet in [37] and will be clarified below. 

As suggested by the transport features discussed 
above, it may not be necessary to use a very long (semi- 
infinite) section of semiconducting material in two doped 
contacts to obtain a strong NDC effect in the P-N junc- 
tions. A finite semiconducting GNM section is expected 
to be enough while limiting the influence of disorder ef- 
fects. To examine this idea, we consider the P-N junc- 




n-doped transition p-doped 




Bias voltage (V) 

FIG. 6: (a) Schematic of simulated graphene p-n junctions 
with finite GNM sections inserted in the two doped contacts, 
(b) shows the obtained I- V characteristics with the effects of 
different Lqnm and disorder strengths. The zigzag holes with 
Q — 15 are used while Ef = and Uo = 0.5 eV. 



tions schematized in Fig. 6(a) and plot their I-V char- 
acteristics in Fig. 6(b) for different lengths of GNM sec- 
tions and for both two cases of with and without disor- 
der. The devices considered here are similar to that in- 
vestigated sec. 3.3 of ref. [37] where a normal graphene 
section is inserted in the transition region, but the very 
long (quasi-infinite) GNM sections in the contacts are re- 
placed by the finite ones. Indeed, without the disorder, it 
is shown that (i) when inserting finite GNM sections into 
the doped contacts, a strong NDC behavior (compared 
to the normal graphene p-n junction) is clearly observed; 
(ii) the PVR of NDC behavior increases and saturates 
when Lqnm is long enough, e.g., Lgnm > 26 nm for 
the devices considered here. It is worth noting that a 
strong NDC behavior with PVR of ~ 200 is observed 
in Fig. 6. However, similarly to the features observed 
above, we find that the current is significantly affected 
by the edge disorder of holes, i.e., though a high PVR is 
still possible, the current peak is strongly reduced when 
increasing the GNM length. In particular, the PVR is 
about 4, 17 and 70 while the observed current peak is 
92, 41, and 3.2 fiA/fim for L GNM = 12.8, 25.6 and 51.1 
nm, respectively. On this basis, to obtain a high PVR 
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and, simultaneously, to avoid the strong suppression of 
the current due to the disorder, it suggests to use (ap- 
propriately) finite GNM sections in two doped contacts. 
For instance, Lqnm — 25^ 52 nm may be a good choice 
for the devices simulated in Fig. 6. 

IV. CONCLUSION 

In this work, using atomistic quantum simulation 
within a tight-binding model solved by the NEGF tech- 
nique, we have investigated the properties of energy gap 
of GNM lattices and considered the optimum graphene 
devices for electronic applications using finite GNM sec- 
tions including the edge disorder of holes. We found that 
the dependence of the GNM energy gap E g on the lattice 
symmetry (i.e., the lattice orientation and the hole shape) 
is significantly suppressed and the quantum-confinement 
effects in the nanoribbon crossing network of GNMs be- 
come dominant in the presence of the disorder. When a 
strong disorder occurs, the dependence of E g on the neck 
width is fitted well with the scaling rule (E g = a/W nw ) 
observed in experiments [26] . which is an important re- 
sult. Regarding the operation of graphene devices, our 



study focusses on single potential barrier structures and 
p-n junctions. It demonstrates that to obtain good per- 
formance (i.e., high ON/OFF current ratio, strong cur- 
rent saturation and NDC behaviors), the use of finite 
GNM sections in the devices can be a very good op- 
tion. More important, since they are unavoidable, the 
proposed configurations are shown to be able to avoid 
the strong disorder effects which tends to limit the out- 
put current through the devices. Though more accurate 
calculation models (e.g., self consistent simulations with 
solving the Poisson equation) are needed for a more quan- 
titative evaluation of device performance, our work pro- 
vides a good guidance for further investigations of GNM- 
based devices. 
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